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Abstract 

In a recent paper [J.Math.Phys. 42, 2236-2265 (2001)], we discussed differential 
operators within a quaternionic formulation of quantum mechanics. In particular, we 
proposed a practical method to solve quaternionic and complex linear second order dif- 
ferential equations with constant coefficients. In this paper, we extend our discussion 
to real linear quaternionic differential equations. The method of resolution is based on 
the Jordan canonical form of quaternionic matrices associoatcd to real linear differential 
operators. 

MSC. 15A18 - 15A21 - 15A30 - 15A90 - 47D25 - 47E05. 

I. Introduction 

Before going into the details of the discussion of R linear quaternionic differential oper- 
ators, we briefly recall the technique used to solve EI and C linear differential equations 
with constant quaternionic coefficients and show the difficulties in extending the method 
of resolution to the R-linear case. We do not attempt a formal discussion of quaternionic 
theory of differential equations. Instead, we take an operational and intuitive approach. 
For the convenience of the reader and to make our exposition as self-contained as possible, 
we follow the mathematical notation and terminology used in our previous paper [1] . In 
particular, the operators 

L^ = (1, Li, Lj, L k ) and = (1, R it Rj, R k ) , \i = 0,1,2,3, 

will denote the left and right action on quaternionic functions of real variable, ip(x), of 
the imaginary units i,j and k. To shorten notation, we shall use the upper-script X = 
{R, C,H} to indicate the R, C, and H linearity (from the right) of n-order quaternionic 
differential operators, 2?*, and n-dimensional quaternionic matrices, M*. 

The general solution of n-order C linear homogeneous ordinary differential equations 



(1) 



where 



v- 



p=0 



has the form 



<p(x) = (fim(x) c„ 



(2) 



where {(^(a;) , . . . , tp 2 n{x)} represent 2n quatcrnionic particular solutions linearly inde- 
pendent over C(l, i) and c m are complex constants determined by the initial values of 
the function >p{,r.) and its derivatives 



ax 

The solution of Eq. (1) is given by 



d ip 



(3) 



(x) = {cx P [M„ c (* - x )] } ^ = £ exp [j„ c {x x )] tp^ , 



where 



P =i 
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and represents the Jordan canonical form of M^. Note that = + where 
is diagonal and N'^ is nilpotent. The Jordan matrix can be determined by solving 
the right complex eigenvalue problem [2-4] for M^. 

It is worth pointing out that from Eq. (4) we can recover the form of the particular 
solutions <p m (x). For example, in the case of a null nilpotent matrix, the general solution 
(2) can be rewritten in terms of left acting quaternionic coefficients (u p and v p ), complex 
exponentials (exp [z p x] and exp [w p x]), and right acting complex constants (c p and c p ) 
determined by the initial conditions. Explicitly, we find 



ip(x) = ^ {u p exp [ z p x ] c p + v p exp [ w p x } c p } , 



(6) 



P =i 



where the complex coefficients {z 1 , w 1 , . . . , z n , w n } represent the right eigenvalues of 
the quaternionic matrix . In the case of equal complex eigenvalues and equal quater- 
nionic coefficients the remaining particular solutions arc determined by using the nilpo- 
tent matrix N%. 

For H linear quaternionic matrices, M^, it is possible to show that v p = u p j and 
VJ V = z p . Consequently, for H linear quatcrnionic differential operators, the general 
solution (6) reduces to 

n n 

<P( X ) = ^2 u p CX P [ z p x ] ( c p + jc P ) = ^2 exp [ q p x } hp , (7) 

p— 1 p— 1 

where q p — u v z v u~ x . The initial conditions (3) shall fix the n quaternionic constants h p . 



Due to the R linearity, the general solution of n-order homogeneous ordinary differ- 
ential equations with quaternionic constant coefficients which appear on the left and on 
the right has the form 



4n 

=^2 ( Ps( X ) r s > ( 8 ) 

s = l 

where {tpi(x) , ■ ■ ■ , (f4, n (x)} represent 4n quaternionic particular solutions, linearly inde- 
pendent over R, and r a are real constants fixed by Eqs. (3). The question still unanswered 
is how to determine the particular solutions tp 3 (x). The natural choice of left acting 
quaternionic coefficients and real exponentials 

<p 3 (x) = u a exp [X a x] , 

does not represent a satisfactory answer. In fact, such particular solutions are at most 
valid for the real part of the eigenvalue spectrum of the matrix M*. 

II. Real linear quaternionic differential equations 

Let us consider the second order R linear quaternionic differential equation 

=0, (9) 

where 

V m - - o (1) — - a (0) = - V a (1) L R — - V a (()) L R 

2 ~ dx 2 dx % - dx 2 ^ ^ ^ dx 2^ <V ■ 



By introducing the R linear quaternionic matrix 



M 2 K = 
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(0) 
T8 


a (1) 



and the quaternionic column vector 



$ = 



we can rewrite Eq. (9) in matrix form 

±<£(x) = M^(x) . (10) 

The real matrix counterpart of the quaternionic operator M", from now on denoted by 
Mg[R], has an 8-dimcnsional eigenvalue spectrum characterized by real numbers and/or 
complex conjugate pairs (the translation tables are given in the appendix). Let Jg[C] be 
the complex Jordan form of M 8 [R] , 

J 8 [C] =R 2m ®Z n ® Z n , m + n = 4, (11) 

where R 2m and Z n represent the matrix Jordan blocks containing, respectively, the real 
and complex eigenvalues of A/ 8 [R] . By using an appropriate similarity matrix 
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<$2m X 2m [R] 


c(4) 
lJ 2mXn 
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£>2m X n [C] 


S„x2m M 


c(5) 
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sj?L [c] 




c(6) 
^nXn 
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we can rewrite M 8 [M] as product of three complex matrices, that is S S [C] J 8 [C]S'^" 1 [C]. 
The problem is now represented by the impossibility to translate the single terms of the 
previous matrix product by R linear quatcrnionic matrices. To overcome this difficulty 
we introduce the complex matrix 



WJC1 = 1, 



1 

V2 



(13) 



An easy algebraic calculation shows that 



t 8 [k] Mmr-'i 



where 
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and 



J 8 [R}=W S - 1 [C}J S 



i?2 



Re[Z n ] 
Im[Z n ] 



-Im[Z„] 
Re[Z n ] 



(14) 
(15) 



The real Jordan canonical form J 8 [R] can be decomposed into the sum of three commuting 
real matrices, that is the diagonal matrix 



D S [R] = Diag |Aj , ... , A m , Ai , ... , A m , Re [zj] 
the anti-symmetric matrix 



Re[z„] , Re[zi] , Re[z„]| 



= 2 







[(! 





Diag {Im 



2l 



Im [z n ]} 



and the nilpotent matrix iV g [R] (a lower triangular matrix whose only nonzero elements 
are ones which appear in some of the sub-diagonal positions). The real Jordan form 
J S [R] = D 8 [R] + A 8 [R] + N a [R] and the similarity transformation T 8 [R] can be now 
translated in their real linear quaternionic counterparts Jf = D^+A^ + Nf and Tf. The 
matrix solution of Eq.(lO) then reads 

$(x) = exp[M 2 R (x - x )]$(x ) 

= If exp[J 2 R (x - x )] (T*)" 1 $(x ) 

= T R cxp[D K (a; - x )] cx P L4 M (x - x )] exp[7V K (x - x )] (T 2 R ) _1 $(x ) . (16) 

In the case of a null nilpotent matrix (explicit examples are given in subsection A), it 
can be verified by direct calculations that the general solution of real linear quaternionic 
differential equations with constant coefficients is 

m 

(p(x) = |n r exp[A r x] a r + u r exp[A r x] /3,.| + 

{v p cos [b p x] — v p sin [b p x] } exp [a p x] 7 P + 



n 

E 

P =i 



{-Dp cos [6 P as] + v p sin [6 P x]} exp [a p x] S p , m + n = 4 , (17) 

P =i 

where a p = Re[z p ], fe p = Im [z p ] and u r ,u r ,v pi v p G H. The 8 real constants a r , (3 r , j p , 
and 5 p are fixed by the quaternionic initial conditions ip(x ) = ip and (x ) = y>i. The 



important point to note here is that the particular solutions corresponding to the complex 
eigenvalue z = a + ib arc given by a quaternionic combination of cos [bx] exp [ax] and 
sin [b x] exp [a x] , namely 

{v cos [bx] — v sin [b x] } exp [a x] and {v cos [bx] + v sin [b x] } exp [a x] . 

For C linear quaternionic differential operators (an explicit example is given in subsec- 
tion B), the general solution (17) reduces to 



(f(x) = u r exp[A r x] c r + v p exp [z p x] d p , c r , d p € C . 



(18) 



A. Null nilpotent matrix 

Let us solve the second order R linear quaternionic differential equation 



dx 



LiRj — LjRi 

dx 



tp(x) = , 



(19) 



with initial conditions ip(0) = j and jjf (0) = k. The matrix operator corresponding to 
this equation is 

1 

LjRi LiRj 



Mi 



(20) 



By using the matrix T 2 

{T^n = (Li + L k Rj — LiR k — L k Ri) — 



1 4 ^ (1 — LjRj + LjRi — R k ) , 



{7?}i2 = ^| (1 + LA - LjRi -L k )-^ (Li + R Z + L } + L k R t ) 



{^?}2i — 2 (LiRt + L k R k — LjR k + Ri — Lj 
{?2 l }22 = (LjR k — L k Rj — LiR k + Rj) , 



Ri 



LiRj) 



and its inverse 

{(T^)^ 1 } 11 = -^j= (1 — LjRj + L k Rj — Li — L k Ri — L,i? fe + LjRi + R k ) , 

|(y*) 1 } 12 = 5 ^q 5 (LiRi + L k R k + LiRj — L k ) + 5 40 - (Lj — Rj — LjR k — Ri) , 
{(T*) _1 } 2 i = (1 + LiRi — LjRi + L k ) , 

{(T 2 K ) _1 } 2 2 = (LjRj - L k R k + LiRj + R k ) + (L k R s - L 3 R k + L 4 R k + Rj) , 

we can rewrite in terms of 



d r = i ( LA + V5L k R k \ a r = ^ 







LjRj 





Ri 



From Eq. (16) we obtain the following solution 



{T 2 K }n exp 
{T 2 K } 12 exp 
I sinh 



LjBi + V5L k R k 



LiRi + ^JZRi 



[m)- i } 1 ij + mr i }i 2 k] + 

[{(r 2 R )- 1 } 21 . ? + {(T 2 R )- 1 } 22 k] 



2\/5 V 



V5 



— y/Z cosh 



sinh 



V5 



2 



V5 
2 

V3 „ 



1 + k 
V3 



V3 



V3 



v3 

2 



| exp [x/2] + 
) } exp [-x/2] . (21) 



In the case of a null nilpotent matrix, the general solution of real linear quaternionic 
differential equations can be written in terms of the eigenvalue spectrum of M 8 [R] , real 
matrix counterpart of M*. In this particular case, the eigenvalue spectrum is 



I 1 
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2 ' 2 
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l-iV3 l + iV3 \ 
2 ' 2 J 



Consequently, Eq. (17) becomes 

1+V5 



ip(x) = Ui exp 
|ui cos 
|i> 2 cos 



VS. 



ai + u 2 exp 



i-VE 



v 2 sin 



■ x 



} 



x a 2 + Ui exp 
cos 

exp [-x/2 i v 2 cos 



| exp [x/2] 71 + |«i 
/2] 72 + {t 



1+V5 



V3 . 



/Si + w 2 exp 



02 



VS. 



+- Vi sin 
v 2 sin 



v3 



I exp [x/2] S 1 
I exp [-x/2] 5 2 



By direct calculations, we can determine the quaternionic coefficients u 12 , u 1 2 , «i. 2 and 
v 1 2 ■ We find 
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61 
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Vi 
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+ k , 
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= i 


+" j , 


"2 
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The particular solutions corresponding to the complex eigenvalues of M s [R] are thus given 

by 



cos 



i±iVa = 
2 ' Zl 



(1 + jfe) exp [x/2] 



V3, 



and 



'-2 



l-i\/3 



z 2 : (« + j) exp [-x/2] 



sin 



2 



2 



The initial conditions y(0) = j and -j^(0) = fc make the solution completely determined, 
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0- 3 
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Let us now consider the following second order R linear quaternionic differential equa- 



tion 



d d 
— 2 ~ (LiRi + LjRj) — — \- Rj., — 



dx" v ' ' J J/ dx 
The eigenvalue spectrum of the real matrix counterpart of 

1 



tp(x) = 



M 2 R 



(22) 



(23) 



is 



{2, -2, 0, 0, 1 + *, -(1-*), 1-*, -(! + *)> ■ 

In this case cosx exp[x] and v 2 cosx exp[— x] does not represent (non trivial) particular 
solutions. In fact, in order to satisfy Eq. (22) we have to impose the following constraints 
on the quaternionic coefficients ?j 1j2 

[L^ + LjRj - 2] v 1 = => V! = i (5 + j 7 , 

[L t R t + L,R 3 + L k - R k ] v 1 = => v x = , 

[L A + Lj^ + 2] u 2 = => i>2 = a + j 5 , 

[L,R, + L,R 3 - L k + R k ] v 2 = => v 2 = . 

The particular solutions corresponding to the complex eigenvalues Z\ = 1 + i and z 2 = 
— (1 + i) are 

j cos [x] — i sin [x] 
Zi , Zi : exp [x] , 
i cos [x] + j sin [x] 



and 



z 2 , z 2 : exp [—x\ 



j cos [a; J + i sm |xj 



i cos [a; J — j sm [x\ 



The general solution is given by 

<p{x) = exp [— 2 x] a.\ + fc exp [2 x] a 2 + fli + fc f3 2 + 

{j cos [x] — i sin [x]} exp [x] 7i + {i cos [x] + j sin [x]} exp [x] S 1 + 
{j cos [x] + i sin [x]} exp [— x] 72 + {i cos [x] — j sin [x]} exp [— x] <5 2 



B. Complex linear case 

We now determine the solution of the second order C linear quaternionic differential 
equation 



77 L-i Rj 

ax 



ip{x) = 



The eigenvalue spectrum of the real matrix counterpart of 



1 

LjR t 



(24) 



(25) 



is 



{1, -1, 1, -1, i, i. 
The general solution is then given by 



i} 



ip{x) = Ui exp [x] «i + u 2 exp [— x] a 2 + Ui exp [x] f3 1 + u 2 exp [— x] /3 2 + 
{v 1 cos [x] — v 1 sin [x]} 7x + {v 1 cos [x] + «i sin [x]} S 1 + 
{v 2 cos [x] — v 2 sin [x]} 7 2 + {v 2 cos [x] + v 2 sin [x]} S 2 . (26) 

By direct calculations, we can determine the quaternionic coefficients u lj2 , w 12 , u 12 and 
u 1)2 . We find 

Ui = fc - 1 , u 1 = j -i , 
u 2 = fc - 1 , u 2 = j - i , 
Vi = 1 + k , «i = , 
w 2 = i + j , v 2 = . 

Consequently, Eq. (26) becomes 

y(x) = (k — 1) [exp [x] Ci + exp [— x] c 2 ] + (f + fc) {exp [ix] di + exp [—ix] d 2 } . (27) 



The initial conditions ^(0) = j and (0) = fc fix the following solution 

i/?(x) = {(j — 1) coshx + (fc — f ) sinhx + (i + j) exp [— ix]} /2 . 



(28) 



C. Non null nilpotent matrix 

As a last example we consider the differential equation 



~T~2 - ( L i R j + LjRi) 

ax 



tp(x) = . 



The eigenvalue spectrum of the real matrix counterpart of 



M 2 R 



f 

LiRj + LjRi 



(29) 



(30) 



|0, 0, 0, 0, y/2, -V2, iV2, -n/2 } 



By using an appropriate matrix Tf, we can rewrite M 2 in terms of 



D K = — 



L 3 R 3 



A R — J- 





Li-Ri 



and = \ 



L, 



- L k Rj 




M 2 M = T 2 K (£>!* + AI + iVf ) (T*)" 1 . Notwithstanding its diagonal form, Nf represents a 
nilpotcnt matrix. In fact, an easy algebraic calculation shows that (iV*) 2 = 0. This 
implies in the general solution the explicit presence of the real variable x. Explicitly, we 
have 



ip(x) = «! + ka 2 + (<5i + ka 2 ) x + 



(i-j) { ex P 
(i + j) | cos 



V2x 
V2x 



[3i + exp 
7i + sin 



V2: 



To make this solution determined, let us specify certain supplementary constraints, 
<^(0)=jand $2(0) = *, 



ip(x) — kx + 3-^- cosh 



V2: 



*2x 



By changing the initial conditions, tp(0) = k and ^g-(O) = j, we find 



ip(x) 



fc + ^|sinh 



V2: 



i+j 
2y/2 



sin 



V2: 



(31) 



(32) 



III. Conclusions 

In this paper, we have extended the resolution method of H and C linear ordinary differen- 
tial equations with quaternionic constant coefficients [1] to R linear differential operators. 
By a matrix approach, we have shown that particular solutions of differential equations 
with constant quaternionic coefficients which appear both on the left and on the right 
can be given in terms of the eigenvalues of the matrix M4„[R], representing the real 
counterpart of the quaternionic operators, M%, associated to the R linear quaternionic 
differential equation. In correspondence to the eigenvalue z — a + ib, we have found the 
following particular solutions 



(33) 



{q cos [bx] — p sin [b x}} exp [ax] and {p cos [b x] + q sin [b x] } exp [a x] , g,pC 
When q = p or p = these solutions reduce to 

q cos [b x] exp [a x] and q sin [b x] exp [a x] . 
For C linear differential operators, the particular solutions (33) couple to give 

q exp [(a + ib) x] . 

Our discussion can be viewed as a preliminary step towards a full understanding of the 
role that quaternions could play in analysis, linear algebra and, consequently, in physical 
applications. A complete theory of quaternionic differential operators (as well of the 
quaternionic eigenvalue problem) is at present far from being conclusive and deserves 
further investigations. 
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Appendix. Real translation tables 
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Example. 
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